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A symmetry conserving, non-perturbative treatment based on a variational squeezed 
vacuum state in conjunction with a well-defined class of RPA fluctuations is applied to 
the SU (3) linear-tr-model. It is shown that the Goldstone theorem holds exactly both 
at zero and finite temperature. The approach represents a systematic procedure which 
avoids problems of the Gaussian Functional with the symmetries. 
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Spontaneously broken chiral symmetry and its restoration at finite temper- 
ature and baryon chemical potential is one of the fascinating and still open 
subjects in strong interaction physics. This problem can be addressed within 
effective theories with linear realization of chiral symmetry, such as linear-cr- 
models [Q. Since spontaneous symmetry breaking and its restoration is non- 
perturbative by nature, non-perturbative methods for solving the effective field 
theory are called for. However, in contrast to a perturbation expansion which 
is trivially symmetry conserving, the situation in non-perturbative approaches 
is more subtle. For instance, in a recent treatment of the three-flavor linear- 
cr-model (Nf = 3) at finite temperature B, to leading order in the 1/iV/- 
expansion, symmetry requirements could not be adequately dealt with, even 
though the 1/iVf -expansion is symmetry conserving. As a result the Goldstone- 
bosons were acquiring masses due to thermal fluctuations. 

In the present letter we show that such difficulties can be avoided by using 
the Quasi-particle Random-Phase- Approximation (QRPA) ||. While mixing 
orders in the coupling constants as well as in the l/./V/--counting, this approach 
is nevertheless symmetry conserving, as has been demonstrated in refs. j|, || 
for the two-flavor case, Nf — 2. The starting point for the QRPA is the quasi- 
particle basis, defined through a variational Hartree-Fock-Bogoliubov (HFB) 
mean field. The HFB theory is equivalent to the well-known Gaussian func- 
tional approach (GFA), extensively used for scalar theories || as well as in the 
context of gauge theories 0. The GFA, although very appealing since it tran- 
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scends the semi-classical approximation based on 1 /TV-expansions, suffers from 
the fact that it is not symmetry conserving and hence the Goldstone theorem 
(GT) is violated. This was, in fact, known since the early days of the GFA 
but no answer has been proposed since, without any guidelines, it is practically 
impossible to guess the missing contributions, necessary to preserve the sym- 
metry. As will be demonstrated below, the QRPA provides a precise framework 
for obtaining these contributions. 

To study SU(3) x SU(3) chiral symmetry non-perturbatively (both at zero 
and finite temperature) we start from the SU (3)-version of the linear-cr-model 
Lagrangian || 

C = ^Tr(d fi Md f *M + ) - ^-TrMM + + f x (TrMM + ) 2 

+ f 2 Tr (MM+MM+) + g (detM + detM+) - c s - c 8 s 8 , (1) 

where /i, f 2 and g are the three possible couplings, [i is the bare mass and M 
the matrix defined by M — -^Aj (sj + ipj). For the Gell-Mann matrices Xj 
our conventions follow those in ref. In the following, sums over repeated 
indices are assumed. For later purposes we recall the conserved axial charge 

Qba = J d 3 xd abc [s b (x)d°p c {x) - p b (x)d°s c (x)] , (2) 

where d a bc is the symmetric tensor, d abc — jTr (A a {Ab , A c }). 

When dealing with the theory at finite temperature and in chemical equilibrium 

the relevant quantity is the grand canonical potential 

n = (H)- TS, (3) 

where T is the temperature, S is the entropy of a gas of interacting bosons and 
(H) the thermal average of the Hamiltonian on the grand canonical ensemble. 
To define the Hamiltonian a second-quantized formulation will be used. There- 
fore the scalar (s^ and pseudo-scalar (pi) fields are represented respectively by 
their creation and annihilation operators bi,bf and ai,af. In a first step, a 
canonical thermal Bogoliubov transformation Q is performed by introducing a 
set of quasi-particle creation and annihilation operators through the following 
rotation 

< T («) - U^^)[uj(q)at( q )-vj( q ) aj {- q )], 

Pt T (l) = U^id^^j^btW-yjiqM^-wjSiq)] . (4) 

where uj(q), vf(q), xf(q) and yf(q) are even functions in the three- momenta, 
and wf denote c-numbcrs defined as 
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The presence of the condensates renders the second Bogoliubov transformation 
in (|4|) inhomogeneous. The orthogonal matrix U decouples those modes which 
are not mass eigenstates of the Hamiltonian, and is given by 

u %i\Q*) = 5 pi s u + ( s i8$ P 8 + S i0 Spo) cosflj + (S iS S p0 - 5 i0 S pS ) sin 6%. (5) 

Here i runs from 1 to 7. Greek symbols will be used to identify the states 
which are mass eigenstates of the quadratic part of the Hamiltonian, and Latin 
symbols for those which are not. However, both letters will designate the same 
set of numbers ranging from 0,..8. The symbol $ denotes both the scalar 
and pseudo-scalar fields and will be used whenever there is no ambiguity. The 
condition, 

W P( g )-,P(g)=*P(g)- y P( g ) = l, 

which renders the transformation in (|j) canonical holds also here. 

In thermal equilibrium, the distribution of maximum entropy is the one which 

minimizes Q with the entropy for a gas of different species of bosons given by 

S = k B ^2[(l + fu)ln(l + fu)-fJnf v ] , (6) 

where /„ are the equilibrium Bose-distribution functions and the sum over v 
includes the number of species as well as the three-momentum q. 
The expectation value of the Hamiltonian in the grand canonical ensemble can 
be easily worked out by using the Bloch-Dominicis theorem [jlO| . This leads to 

(H) = J ^u Jq [{l + 2fl{ q )){u 2 M)+vUq)) + (l + V!{q)){x 2 M+vK<l))] 

1 

+ c i£i H _ Gijkiiij£,k ~ -^^ijkl^j^k^l ~ (^ijkl^k^l + 3Gijk£k) ^ij 



where Fijki, Fij,kU Gijk are tensors which have also been given in [ pT[ 

Fijkl — — (3/l<5 m O + /2) [dijmdklm + dumdjkm + dikmdjlm] 
J~ ij,kl — 2 v^fl^mO ~t~ f^) [dijmdklm ~t~ filmfjkm fikmfjlm] 



(7) 



zfijk — "^3 
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2 {$i0 + $jo + Sko) + 2 SjoSko 



(8) 



To fix the functions Ui(q), Xi(q), 9 Sl 9 p and Wi we make use of the Ritz varia- 
tional principle. Minimizing the grand potential f2, the finite-temperature BCS- 
equations follow straightforwardly. The quasi-particle energies for the scalars 
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in terms of the scalar and pseudo-scalar gap parameters 

3 



(9) 



Fijklikil + -jQijk^i + Fijkll{ s)kl + J r ij,kl I (p) k 

rp rp 3 T T T 

^ijM^k^l ~ T^Sijk^k + -Fijkll( p ) kl + ^ij,kl^(s) kl 



(10) 



The one-point functions are denoted by = (ui), while ^^v. are thermal 
averages for the bilinears of pseudo-scalar and scalar fields. At the minimum, 
these can be expressed explicitly in terms of the U -matrices as 

= / dx = C/ (*) 1P C/ (*) J /(*) P > (H) 

where 1^ are thermal tadpoles given in terms of the self-consistent quasi- 



particle energies and the Bose factors (q) by 
d 3 q l + 2/ ( T $) >) 



(2tt)3 2£ ( ^ )p (q) ' 



wit/i /($) p (g) 



exp 



(12) 

It is important to notice that this self-consistency is not forced by hand but 
rather follows directly from the Ritz principle, as was shown in g g. The 
mixing angles are fixed through the equation 



A (*)ij i^oSjo - SisSjs] sin26»$ T = Aj^ y [5 i0 S j8 + 5 iS S j0 } cos20$ T . (13) 



Finally, the minimization with respect to Wi leads to an additional condition 
which will determine the strange- and non-strange condensates, related to £o = 
(00) and £ 8 = (a s ) via 

+ {&mi(j ~ 3&fcl) I( P ) kt for i = 0or8. (14) 

The BCS equations (|[ |l3[ |l4|) form twelve coupled, self-consistent equations. 
At zero temperature, they completely define the variational ground state of the 
theory \^>) as a squeezed state 1 . In the context of functional field theory ap- 
proaches this ground state is equivalent to the Gaussian functional which has 

1 The squeezed vacuum is implicitly defined by the relations : ap(g)\^f) = (3 p (q)\'i/) = 0, for 
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been actively studied in the context of the triviality of A<I> 4 theories |T^, 
According to these studies, a spontaneously broken phase in the vacuum of 
the renormalized 0(iV)-vector model can be secured. This is a very interesting 
feature, since it is ruled out in the 1/A^-expansion, where the vacuum of the 
(renormalized) theory lies in the same group of invariance as the Lagrangian. 
An additional complication arises, however, in case of broken continuous sym- 
metries. The squeezed vacuum is not compatible with the Goldstone theorem 
which calls for the appearance of as many Goldstone modes as there are broken 
generators. This is clearly not the case here, as was shown in many instances 
(see [|], [l4], |l5|). We will see below how this can be corrected. 
It is well known that the self-consistent HFB mean field is a 'mixing-order' ap- 
proximation (see for instance ||) in the sense that there is no ordering in terms 
of the coupling constants nor in the number of flavors (colors or charges in other 
situations). At first sight this leads to difficulties, as there are apparently no 
guidelines for classifying the relevant missing contributions for preserving the 
symmetries. A powerful approach, used for the first time in the context of the 
SU(2) linear-c-model in ref. Q, has proven successful in this regard. We will 
see that it is capable of handling the much more complicated 5t/(3)-case as well. 
The basic idea for restoring the symmetry (broken at the HFB-level) is to use 
the symmetry generator for defining a general RPA excitation operator. As is 
well known in many-body physics, RPA fluctuations, built on a self-consistent 
mean-field, restore broken symmetries. If the latter are continuous this is real- 
ized through the appearance of spurious excitations of zero energy in the RPA 
solutions. We briefly sketch how this can be achieved. First one recalls the RPA 
equations which, for the finite temperature, read |l6| 

([SQ U , [H,Q+]])=m v ([6Q v ,Q+]) , (15) 

where the expectation values are taken in the grand canonical ensemble and 
m„ denotes the excitation energy (mass) of a given state. The excitation op- 
erator Q+ is chosen so as to generate pseudo-scalar excitations corresponding 
to the 7r, K , K, rj and rf mesons. The RPA equations can be derived from a 
variational principle in a restricted Fock space. A further approximation is em- 
ployed, the so-called quasi-boson approximation ( [Q v , Qt\ ) = 8 Vfl , which allows 

any state p. It is obtained by successive unitary transformations to the trivial vacuum, such 
that: I^P) = UmxUshUsq Hi- The first two U sq and U a h are respectively a squeezing and a shift 
unitary transformation 17[, while U m x is a mixing unitary transformation |18| . Since these 
three can be expressed as exponentials of, at most, bilinear forms in the creation and annihila- 
tion operators, the vacuum \ ^) can therefore be put in a more compact [JL9|, but less appealing 
form which, in the case of a SU (3)-invariant vacuum, as treated at the end of this paper, reads 

E 9 5$&4to>4(-«> + £jrh b l^ b l(-i) + |o>. 



\\S/) = A/"exp 
Af is here an irrelevant normalization factor. 
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for a linearization without destroying the symmetry conserving properties of the 
thermal RPA. The RPA equations ( |l5| ) will generate spurious solutions if the 
operator contains the symmetry generator specified by Q^ a in @- In the 
chiral limit (cq = cs = 0) the Ihs of ( |l5[ ) will then be zero since the Hamiltonian 
commutes with the symmetry operator. If the norm on the rhs is well behaved 
then this forces m v to be zero. Normalizability is guaranteed, however, by the 
use of the HFB quasi-particle basis. Hence m v vanishes in the chiral limit at all 
temperatures as long as the symmetry is spontaneously broken. In the following 
this will be shown explicitly. 
To proceed we take Q+ as 

Q+ = ^Ka+(0)-^ flp (0)] 

+ M p (q){3+(q)a+(-q) - V^ p (q)^(q)a p (-q)] 

i 

+ ^^[^(^feKd)- ^(^(!)a+(4)] , (16) 
i 

where U, V, W, X, Y and Z are general amplitudes which will be fixed once 
the RPA eigenvalue problem solved. Expressing Q$ a (||) in second-quantized 
form it is evident that Q„ contains the same excitations as Q^ a 2 . There is 
an interesting aspect which is apparent from the ansatz ([l6]) of the excitation 
operator. The latter, in addition to the RPA bilinear forms which as usual 
induce quasi-particle scattering processes, contains linear forms which indicate 
that the RPA scattering equation is in fact coupled to a Dyson equation. Note 
also the (crucial) presence of the tensor d vpp which will act as a projector on 
the relevant RPA scattering-states, necessary for fulfilling the GT. It is clear 
that one has a non-trivial criterion for selecting the relevant fluctuations which 
are difficult to guess without having a systematic approach. 
For the purpose at hand, the RPA equations can be condensed as 

^ 1 ^up, an dvoT] Ilq.^ = , (17) 

alj 

where FI^ is some eigenvector which does not need to be specified any further. 
The other quantities are defined as follows 

U (phr, U (s)i„ E ^Pn( m 

a Starting from Q^ a in (H), one first expresses the fields in second-quantized form using 
a + , a, b + , b. Then one inserts the inverse Bogoliubov-transformation (wh. This leads to an 
expression for Q$ a similar to the one for in jl6[ ). However, in this case the amplitudes U, 
V, W, X, Y and Z will have well-defined forms while in the general case of they remain 
to be fixed by the eigenvalue problem (Jl 5[) ■ 



M, 



HP, an 



8F. 



ij,kl 



R Tlk R Tjl 



,T2 



?T1 
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where S SpPij (m^) are scalar pseudo-scalar thermal RPA loops given by 



^p^v ) - J J^y 



£l(q)+£l(q) l + fZ(q) + fZ(q) 
2££ (q)£l (q) m™ ~ (££ (q) + & (q)Y 



£^q)-£iM) fj{q)~UM) 
2£T(q)ET {q) m T2 _ {S T {q) _ S T (q) y 



(19) 



These RPA equations, build on the self-consistent HFB basis, are valid for the 
whole nonet of pseudo-scalars. In the chiral limit, the octet of pseudo-scalars 
ought to come at zero RPA frequency while the singlet is not a Goldstone mode 
as long as the anomaly term is present in the Lagrangian (g ^ 0) . This is indeed 
reflected in the RPA equations since the ninth axial charge does not commute 
with the Hamiltonian as long as g =/= 0. 

To verify that the symmetry is recovered at the RPA level a decisive test is to 
see whether the octet of pseudo-scalars has zero mass when the symmetry of the 
vacuum is broken from the SU(3) x SU(3) group down to SU(3) which happens 
when the condensate (erg) vanishes. The c% explicit breaking is consequently 
put to zero (cs = 0). Close inspection shows that in this situation, as in the 
perturbative case, the mixing angles 6 S and 9 p vanish 3 . Therefore, at the 
HFB level, eqs. (0) generate mass eigenstates of the Hamiltonian in this limit. 
Accordingly, we will not distinguish any longer between states with Greek or 
Latin labelling. Finally, from the HFB equations, one can infer the following 
solution which is consistent with the octet degeneracy: 

£(p) k = £n I £(s) k — £a ' f or k =1, ..8 , 

E l)o = £? s)0 =£j ^£j ■ (20) 

This solution is, of course, compatible with a finite condensate (ctq) and an 
overall mass splitting between the scalars and pseudo-scalars. It is obtained for 
a given choice of the parameters in the Lagrangian. Assuming this solution, 
the bare mass /z can be fixed accordingly. Injecting the relations (^0|) into the 
RPA equations ( p7|jT8| ) one can derive the RPA eigenvalues. After lengthy but 
straightforward algebra which consists essentially in solving the characteristic 

3 Explicitly one has .F08OO = ^"0800 = 0080 = 0. Furthermore Fo8klI(*s>) kl and ^"o8,fc!-f(*) fci 
are both proportional to d% k il^ . On the other hand the tensor alg k [ enforces k = I ^= 
0. Assuming octet degeneracy ( I( p ) kk = Itt, I(s) kk = It f° r k = 1,8) and the fact that 

X/fc— 1 d&kk = 0, dsoo = 0, one finally has a self-consistent solution with A^*'<5io<5;8 = 
and hence 9$ = 0. This solution is consistent with the octet degeneracy of the scalars and 
pseudo-scalars and is reminiscent of the perturbative case. 
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equation of the RPA eigenvalue problem, one can read off the RPA frequencies. 
Denoting the masses of the octet of Goldstone bosons generically by one 
finally has 



T2\ 



with the following definitions 

s T K 2 ) 



1 - |/ 2 [ 2 C« 2 ) - K,(ml 2 ) - Co {ml 2 ) 



32/|S a7r (m^ 2 ) 
5u 2 _t . T 2 n 



T 

Kp (m 



T 2x 



2(u + w)" — i 



T 2\ 



8/2 



-w 2 - (u + w) 2 
9 

^v 2 -{u + vf 



2(u + u) 2 _t 



K 2 ) 



s;(mI 2 )CK' 2 ) 



T 2\ 



^/ 2 (,-») 2 4( m D4( m r) 

64/| 



(« - wfK,{rnl 2 )K P A< *)K A< 2 ) 



T 2 



T 2\ 



(21) 



2V2q - 



,K 2 ) 



,T 2a 



1) = -3V2Q 



Pij K 2 ) l + S/rE^K^) 



-3V2<? - 4V6/1 C T » 
(22) 



As a last step, we show that this expression supports a zero-energy solution. 
We briefly sketch the proof here. Further details will be deferred to j2Cj. 
As long as the symmetry is spontaneously broken, equations (||) for £j, £p g , 
£j, £j and ( [l4|) for ^ allow for a solution with finite £j . Assuming this, one 
can substitute the bare mass /x in (§) by its value from (fj). When taking into 
account the relations in (20) one can express the quasi-particle energies as a 



function of the thermal RPA loops with zero argument: £^(0), E^ po (0) and 
Sj o7r (0). These uniquely specify the solution since the remaining possible RPA 
loop Ej opo (0) can be expressed as a function of these three. Using the identity 
for regularized integrals 



L (p)i 



?T 2 



(0), 



(23) 
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which holds for zero as well as finite temperature, one can verify that the ex- 
pression for the quasi-particle energies is exactly given by 



T2 __A^(0)_cp_ 



This clearly indicates that in the chiral limit (co = 0) and for a broken phase 
(Ctf 7^ 0) a spurious solution in eq. with = exists and remains as 

such, protected by the symmetry, until the transition temperature T c is reached. 
Once the transition to the Wigner-Weyl phase occurs, the spurious solution in 
the RPA disappears due to the fact that the corresponding eigenvector is not 



normalizable . Eqs. (21) and ( |24[ ) are the mathematical reflection of the 
GT. They form in fact the Ward identity which links the Goldstone two-point 
function to the condensate 5 



-D- 1 (0) = 



-c 



One can appreciate the non-triviality of the present result which suggests that, 
for a one-point function (condensate) evaluated at the HFB level, the corre- 
sponding Goldstone mode has to be an excitation in the spectrum of the RPA. 
As mentioned earlier, this result is valid for the octet of Goldstone modes while 
the singlet pseudo-scalar remains massive since the corresponding symmetry 
generator does not commute with the Hamiltonian, because of the anomaly. 
Further details as well as a quantitative study of the phase transition will be 
presented elsewhere [ p0| . 

In summary, we have discussed the SU(3) linear-cr-model at finite tempera- 
ture within the HFB-RPA approach. The Goldstone theorem has been shown 
to hold both at zero and finite temperature below the phase transition. Above 
the transition the spurious solution of the RPA equations are no longer nor- 
malizable, therefore no asymptotic zero-energy solution is present in the RPA 
spectrum. An explicit proof of the Goldstone theorem has been given in the 
case of the vacuum being an invariant of the SU(3) symmetry group. It should 
be appreciated that the present approach is transcending the usual expansion 
in l/Nf. It is in fact order mixing but nevertheless symmetry conserving. It 
shares with the Gaussian functional approach its advantages and corrects for its 
shortcomings when it comes to fulfilling the symmetry requirements. We have 

"^Some of the elements of the diagonal matrix-norm given by the rhs of eq. (|lf^) are propor- 
tional to the differences / Sp ((?) — fp v (<?)• These vanish in the Wigner-Weyl phase due to the 
degeneracy of scalar and pseudo-scalar modes. In this situation the eigenvalue problem in 
(jig ) does not have a solution any longer (see ref. 

^Eq. ( |2l| ) is the pole equation which leads to the pion mass. The residue at the pole which 
receives contributions from the RPA fluctuations, is absorbed in the definition of co. More 
details on Ward identities in the linear-sigma-model could be found in |2l| . 
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also demonstrated that this approach can handle successfully not only models 
with O(N) symmetry W but also those with SU(N) symmetry. 
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